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First principles quasiparticle damping rates in bulk lead 
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First principles calculations of the damping rates (inverse inelastic lifetimes), r^ 1 , of low energy 
quasiparticles in bulk Pb are presented. Damping rates are obtained both for excited electrons and 
holes with energies up to 8 eV on a set of k vectors throughout the Brillouin zone (BZ). Strong 
localization effects in the calculated r _1 are found. Averaged over the BZ inelastic lifetimes versus 
quasiparticle energy are reported as well. In addition, the effect of the spin-orbit induced splitting 
in the band structure on the calculated lifetimes in Pb is investigated. 

PACS numbers: 71.10.Ca.71.15.Mb.75.70.Tj 



I. INTRODUCTION 

In metals electron-electron inelastic scattering pro- 
cesses give rise to the main contribution to the damp- 
ing rate of excited electrons (and holes) with energies 
>0.5 eV above (below) the Fermi level. For a long time 
the basic knowledge for such kind of processes was based 
on theories developed for a free electron gas (FEG).— ~— 
Recently, when the calculations from the first principles 
became to be computationally feasible, this field has ex- 
perienced profound modifications. Thus the first princi- 
ples calculations^— have shown that the inelastic life- 
time of excited electrons indeed is a result of balance be- 
tween localization, screening and band structure details, 
even in metals whose electronic structure frequently con- 
sidered as being a free-electron-like onei 15 ' 16 In the pre- 
vious theoretical work a non-free-electron like behavior 
of damping rates in Be ; 15 ' 17 the role of the screening of 
the d-electrons in inelastic lifetimes on Cu^^ transient 
excito n 18 ' 19 and full inclusion of exchange-correlation 
(XC) effects 2 ^ as well as renormalization effects on quasi- 
particle (excited electron or hole) lifetimes in noble and 
transition paramagnetic^— and ferromagnetic metals^ 
and compounds^ 5 - - — were studied by means of the GW 
approximation. 28 Using the T-matrix theor y 29 ' 30 the con- 
sequences of spin-flip processes on damping rates in mag- 
netic materials were found to be significant for the spin- 
minority states. 24 ' 31 ' 32 In Ref. [33| good agreement was 
found between the inelastic lifetimes for excited electrons 
and holes in several metals evaluated within the GW ap- 
proximation and the semiempirical scattering-theory ap- 
proach. Role of inclusion of the accurate quasiparticle 
band structure of Cu and Ag in the quasiparticle life- 
times was investigated . 34 ' 35 

At the same time, in heavy elements an additional 
ingredient like spin-orbit (SO) interaction starts to be 



important in the description of the electronic structure. 
A well-known example is corresponding modifications in 
band structures in Bi and Pb. However, to the best of 
our knowledge, up to now the effect of the SO interac- 
tion on quasiparticle lifetimes in real materials taking 
into account its band structure evaluated from first prin- 
ciples was not investigated. In the present paper, for 
the first time, inelastic lifetimes of excited electrons and 
holes are studied in bulk lead by means of first-principle 
calculations, analyzing in detail the band structure as 
well as SO coupling effects. Recently thin films of Pb 
grown on different substrates have obtained great deal 
of attention. The questions regarding growth, transport, 
magnetic, superconducting properties of these systems 
have been considered. For instance, the confinement ef- 
fects on the superconducting transition temperature 3 -^— 
and quasiparticle decay rate s 41 ' 42 were studied. Recent 
lifetime measurements by two-photon spectroscopy 4 ^ for 
few monolayers Pb deposited on silicon substrate point 
out to the conservation of the bulk-like behavior of decay 
processes in Pb even in the very thin slabs. Therefore 
detailed first principles investigation of the inelastic de- 
cay rates of quasiparticles in bulk lead which can serve 
as a reference for existing and future experiments seems 
useful and timely. 

In the present study the Kohn-Sham equations of 
density-functional theor y 43 ' 44 are solved self-consistently 
within norm-conserving pseudopotential scheme. Sub- 
sequently, linear response theory is used to calculate 
the momentum- and energy-dependent density response 
functions over the Brillouin zone (BZ) from which the 
imaginary part of the quasiparticle self-energy is com- 
puted invoking the GW approximation of many-body 
theory. Comparison of the obtained results with that 
for the FEG model shed light on the role of the screen- 
ing in this metal. In order to study the effect of the 
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SO coupling on the quasiparticle lifetime through its in- 
fluence on the band structure the latter was evaluated 
both including and excluding the SO term in the Hamil- 
tonian. The resent work demonstrates that the inclusion 
of the SO interaction into the band structure calculation 
produces noticeable effect in inelastic lifetimes for quasi- 
particles at very low energies. Nevertheless, the general 
trends in the quasiparticle inelastic lifetime as a function 
of its energy are unaffected by the inclusion of the SO 
splitting. Additionally, comparison of the calculated in- 
elastic lifetimes for the p and d electrons with the same 
energy reveals strong localization effects. 

The paper is organized as follows. In Sec. II, the 
explicit expressions for the quasiparticle decay rate in 
periodic crystals within the GW approximation are pre- 
sented. In Sec. Ill the effect of the SO interaction on the 
electronic band structure of bulk Pb is discussed. The 
calculated results for the damping rates are analyzed in 
Sec. IV along together with a link to the linewidth of the 
quantum well states (QWS) in Pb(lll) thin films. Fi- 
nally, conclusions are given in Sec. V. Atomic units are 
used throughout, i.e., e 2 = ft = m e = 1, unless otherwise 
is stated. 



II. CALCULATION DETAILS 

In the framework of many-body theor y 29 ' 30 for an in- 
homogeneous electron system, the damping rate rf of 
an excited electron (or hole) in the state </>i(r) with en- 
ergy Si is obtained from the knowledge of the imaginary 
part of the quasiparticle self-energy, S(r, r'; e^), as 



-2 Jdr y"dr>*(r)Im^r,r';£. i )<Mr')- 



(1) 



In order to calculate r _1 via Eq. ([T]), the nonlo- 
cal energy-dependent operator S(r, r'; e^) is evaluated by 
means of the GW approximation^ 2 — in which one consid- 
ers only the first order term in a series expansion of the 
self-energy in terms of the screened Coulomb interaction: 

S(r, r'; e t ) = ±- J G(r, r'; £l - s)W(r, r'; e)de, (2) 

where G(r, r';Si — e) stands for the one-particle Green 
function and W(r, r';e) is the dynamically screened 
Coulomb interaction. Replacing the Green function by 
its noninteracting counterpart, the imaginary part of the 
self-energy can be expressed in the following way: 

ImS(r,r'; £i ) = ^^(r')ImW(r,r';a;)^(r), (3) 
/ 

where w = \e% — £/| represents the energy transfer and the 
sum is extended over a complete set of final states <fif(r) 
with energy e f (with constrain Ep ^ £/ ^ £i for electrons 
and E{ ^ Sf ^ Ep for holes), being Ep the Fermi energy. 



The screened Coulomb interaction is written as 

W(r, r'; lo) = v(r — r') + J dri J dr 2 v(r — ri) 

x x(n,r 2 ; uj)v(r 2 - r'), (4) 

where v(r — r') is the bare Coulomb interaction, and 
x(r,r';w) is the density response function for interact- 
ing electrons. 

In the framework of time-dependent density- functional 
theory 4 -»^& and within the random-phase approximation 
(RPA), the density response function x satisfies the in- 
tegral equation 



X( r , r ';u) = x°(r,r';uj) + J dr x J dr 2 x°(r,r 1 ;uj) 

xv(r 1 -r 2 )x{r2,r';oj), (5) 

where x°( r : r '; w ) is the density response function of the 
noninteracting electrons. Choosing the RPA one neglects 
the short-range XC effects in the evaluation of x(r, r'; tS). 
At T=0 K, x°(r,r';w) rea ds 

f£ £j ~ ej> + (oj + u]) 

x^(r)0*,(r)^(r')^(O- (6) 

For tree-dimensional periodic crystals the following 
Fourier expansion for the screened interaction of Eq. ^ 
can be introduced: 

1 BZ 



q g,c 



xu G (q)e G . G ,(q,w), (?) 



where the first sum is extended over the first BZ, G 
and G are reciprocal lattice vectors, f2 is the normal- 
ization volume, v G (q) represents the Fourier coefficients 
of the bare Coulomb interaction, and e G 1 G /(q, w) are the 
Fourier coefficients of the inverse dielectric function, 

e G * G ,(q,w) =<J G ,G' +XG, G '(q,w)v G <(q). (8) 

Introducing Eq. ((7} into Eq. the following expres- 
sion for the damping rate of a quasiparticle in the initial 
state ni k(r) with energy e ni k is obtained: 

v^E^E ^ 

xIm[-e G i G ,(q, W )], (9) 

where u = e Ulk - e n/ k-q, and 
B„ in/ (k,q,G) = J dr^ ik (r) e ^ +G )- r 0„ /k _ q (r). (10) 
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Note that all matrix elements of inverse dielectric func- 
tion q/ (q, uj) enter Eq. © because of the coupling 

between wave vectors q + G and q + G' with G 7^ G' as 
a consequence of the electron density variation in solids. 
These coupling terms represent so-called crystalline local- 
field effects4I 

For evaluation of the density response function 
X°(r,r';w) and the matrices B ni71f (k, q, G) the eigen- 
functions of a Kohn-Sham system without SO interaction 
are expanded in a plane- wave basis, 



Ilk (r) = -^^ Mnk (G)e^+ G )- r . 



(11) 



When the SO interaction is included into the Kohn-Sham 
Hamiltonian, one has to work with spinors as eigenfunc- 
tions. In this case, 



1 



$nk(r) = -=^J]iWk(G; 



= i(q+G)-r 



6 



(12) 



G 



V2O 



for spin up and £4, 



1 

V2 



°) for spin 



where £| 
down. 

In this work the electron-ion interaction was repre- 
sented by norm-conserving nonlocal pseudopotential,— 
and the local-density approximation was chosen 
for the exchange and correlation potentials in the 
Ceperley-Alder forrn^ using the Perdew-Zunger 
parametrization^ Well- converged results for the band 
structure have been obtained with a energy cutoff of 14 
Ry which corresponds to the inclusion of ^180 plane 
waves in the expansion of the Bloch states. 

The calculations of the inverse dielectric matrices 
through Eq. ([5]) were carried out using for evaluation 
of x° m Eq- @ a k sampling over « 32000 vectors in 
the irreducible part of BZ (IBZ) and including 25 energy 
bands. The broadening parameter 6, employed in the 
evaluation of x° as explained e.g. in Refs. [sT| and 
was set to 10 meV. In the expansion of dielectric ma- 
trices 40 plane waves have been considered. The sums 
over reciprocal vectors G and G' in Eq. ^ have been 
extended over 40 vectors as well. 



III. SPIN-ORBIT EFFECTS IN PB BAND 
STRUCTURE 




FIG. 1: (Color online) Calculated band structure of bulk lead, 
with (solid lines) and without (dashed lines) inclusion of the 
spin-orbit coupling. Squares and circles mark regions where 
the SO-induced band splitting is reflected in the density of 
states shown in Fig. [2] The horizontal dotted line represents 
the Fermi level. 



theoretical results 5 ^ and with the experimental data 55 
when the SO term is taken into account. As can be seen 
in Fig. [TJ the inclusion of the SO interaction affects only 
three p energy bands (the only ones that cross the Fermi 
level) mainly around the high-symmetry points. In Ta- 
ble HI a comparison of the energies of the p bands at the 
high-symmetry points with the inclusion and exclusion 
of the SO coupling is presented. One can see that the 
SO interaction produces two main effects in the Pb band 
structure: 

a) At the scalar-relativistic level, at some points and 
some high-symmetry lines of the BZ, the number of states 
in the the energy bands can be more than two (however 
always it must be an even number because of the Kramers 
degeneracy). The SO coupling lifts some of these "ex- 
tra" symmetry degeneracies. The largest SO-splittings 
Ae so are observed for p electron bands at the T, X, L, 
W points and close to K, as well asalong the respective 
symmetry directions. At the high-symmetry points we 
obtain: Ae so = 3.21 eV at the T point, 1.20 eV at X 
and W, and 1.85 eV at L (see Table©. 

b) The SO interaction leads to the avoiding band- 



Bulk Pb lattice has a face-centered-cubic (fee) crystal 
structure. Fig. [T] presents the calculated band structure 
for lead along the high-symmetry directions of the BZ ob- 
tained with the use of the experimental lattice parameter 
a c — 4.95 A. In this figure one set of data (red solid lines) 
corresponds to the calculation with the inclusion of the 
SO interaction in the Kohn-Sham Hamiltonian whereas 
black dashed lines present results with the exclusion of 
SO interaction. As the fee lattice has inversion symme- 
try, due to the Kramers degeneracy 5 ^, each energy band is 
at least double degenerate in spin in both cases. The cal- 
culated band structure is in good agreement with other 



TABLE I: Energies Ei of p-like electronic states at some high- 
symmetry points in the Brillouin zone obtained in scalar- 
relativistic (sc) ab initio calculations and with SO interac- 
tion included computations (so). All energies are in eV with 
respect to the Fermi level. 





r 


X 




W 


L 




K 






sc so 


sc 


so 


sc so 


sc 


so 


sc 


so 


£l 
£2 
£3 


7.14 4.97 
7.14 8.18 
7.14 8.18 


-3.18 - 
1.73 
1.73 


3.30 
1.20 
2.40 


-1.40 -2.09 
-1.40 -0.89 
0.28 0.45 


-4.07 - 
5.23 
5.23 


4.13 
4.30 
6.15 


-2.25 - 
-1.19 - 
2.70 


2.52 
1.04 
2.79 
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FIG. 2: (Color online) Total density of states (DOS) obtained 
in the scalar-relativistic calculation (DOSsc, dashed line) and 
the calculation with inclusion of the spin-orbital (SO) term 
(DOSso, solid line). The main variations in the DOS upon the 
SO inclusion reflected in differential ADOS=DOS S o-DOS S c 
(dashed-dotted line) are related to different band structure 
splitting marked by the same symbols as in Fig. [JJ DOS is in 
arbitrary units and energy is with respect to the Fermi level. 



crossing effect, as the ones between p-like states marked 
by circles in Fig. [TJ 

In Fig. [5] the total density of states (DOS) is plot- 
ted as a function of energy. The most important ef- 
fect of inclusion of the SO interaction is the weight loss 
around -1.4 eV and the increase of the peak at -2.4 eV, 
whose position is also slightly shifted to higher binding 
energies. These variations reflect disappearance of the 
band-crossing points mentioned above. The other effect 
- breaking of the band degeneracies - has no important ef- 
fect on the DOS, because it does not flatten significantly 
the band dispersion. Nevertheless, the SO-induced split- 
ting around X and L points affects the DOS at ^1.6 eV 
and ^5 eV, respectively (dips in ADOS in Fig. f2]markcd 
by squares). 

In the band structure calculations the Hamiltonian in- 
cluding the SO term was solved fully self consistently. 
However, in the lifetime calculations the spinor represen- 
tation, Eq. (fT2"j). for the relativistic wavefunctions is not 
used. First, including the full spinor structure of the rela- 
tivistic wavefunctions dramatically increases the compu- 
tation time. Secondly, it is well known that in simple sys- 
tems with strong spin-orbit coupling, the SO-splittings of 
the one-electron energies are important. Nevertheless the 
total electronic density and furthermore the partial con- 
tribution of each one-electron state to the total density 
show negligible variations upon switching the SO inter- 
action (see, e.g., Ref. [H). Hence, we suggest that one 
could mostly retain the SO effects on r~ k just through 
the inclusion in Eqs. © and (J9j) of the one-electron ener- 
gies calculated including the SO term in the Hamiltonian, 
while using the scalar-relativistic wavefunctions Eq. 
in Eqs. ^ and (fTU|) . As a result, in the SO case the 



energy-loss function Im[— e G G ,(q, w)] entering Eq. © 

for t"^ differs from that used in scalar-relativistic calcu- 
lations. In a future work, lifetime calculations including 
the spinor representation of the relativistic wavefunctions 
will show if this is a meaningful approximation for includ- 
ing the SO effects on the calculated r„ 
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IV. DAMPING RATES 

In Fig. |3] we show the calculated damping rate r„k = 
r nk °f n °t electrons and holes for all the values of band 
index n and wave vectors k used in the calculations, both 
using the scalar-relativistic energy spectrum (black dots) 
and including the SO coupling in the band structure (red 
dots) . In addition, the blue dots represent the results ob- 
tained using the Lindhard dielectric functional on calcu- 
lating the screened interaction through Eq. ([7]) together 
with the ab initio scalar-relativistic wave functions and 
one-particle energies in Eqs. (JlJ and (J3j) . The orange line 
is the result of a FEG model, i.e. when the band struc- 
ture is described by a band with parabolic dispersion, the 
wave functions are represented by plane-waves, and the 
Lindhard dielectric function is used for the screening. In 
the inset of Fig. [31 the green curve represents the Quinn- 
Ferrel result for the damping rate for very low energy 
quasiparticles: 



.5/2 



'QF 



(e - E F f 
263 



eV- 2 fs- 



(13) 



where r s is the valence charge density parameter of the 
system. In the following we shall refer to these five kinds 
of lifetime calculations as SC, SO, LDF-SC, FEG, and 
QF ones, respectively. 

There are three main features in the distribution of the 
damping rates in Fig. [3J First, for -8 eV < e < -6 eV 
the damping rates corresponding to the s-like holes are 
presented. On the opposite side of the studied energy 
range, the values of T n k for excited electrons in the d- 
like bands appear for energies in the 5.5 eV < e < 8 
eV interval. However the main feature of Fig. [3J is the 
parabolic-like distribution of damping rate for the p-like 
holes and electrons in bands crossing the Fermi level. The 
linewidth data for the p states is separated from that for 
the s-like holes by the energy gap of «2.5 eV presented 
in the Pb band structure. 

It is interesting to compare the observed difference on 
the calculated damping rates for p and d electrons with 
the same energies. As can be seen in Fig. [31 at the same 
energy in the 5.5 eV < e < 8 eV interval quasiparticles 
in the d bands present a damping rate roughly two times 
smaller in comparison with that for the p bands. This is 
a consequence of the different coupling matrices [see Eq. 
(10)] for p and d states thereby signalling about strong 
localization effect in lifetimes for ci-states in Pb. 

For quasiparticles in the p-like states, the role of final 
states in its decay [see Eq. (10)] can be seen in Fig. 
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FIG. 3: (Color) Calculated damping rate for hot electrons and holes in all the bands and k points set over the irreducible 
part of the Brillouin zone included in the calculations. Red dots represent the values calculated including the spin-orbit (SO) 
coupling in the band structure, the black lines represent the scalar-relativistic results and the blue lines are the results obtained 
using the Lindhard dielectric function together with the scalar-relativistic band structure. The orange line shows linewidth 
obtained from a free electron gas model with the Pb valence charge density parameter r s = 2.298 a.u. The green curve in the 
inset represents results for linewidth calculated according to the Quinn-Ferrel expression according to Eq. (|13[) f° r the same 
r a . The black arrow points to the linewidth results for the lowest p band near the BZ center calculated with the inclusion of 
the SO coupling term (see text for more details). 



TABLE II: Linewidth for the lowest p states at k's in vicinity 
of the F point along the TX direction, calculated at the scalar- 
relativistic level (F sc ) and including the SO splitting (r so ) in 
the band structure. All values are in meV. The k vectors are 
labeled with the k index, corresponding to k = 0, 0), 

where a c is the lattice parameter. 



K 





1 


2 


3 


4 


5 


6 


■psc 


1135 


1100 


1010 


895 


776 


663 


559 


j-iSO 


440 


430 


410 


381 


345 


307 


268 



[3] from comparison of the distribution of the LDF-SC 
results (blue dots) with the FEG line. As in both these 
cases the Lindhard dielectric function was employed, the 
difference arises only in the wave functions entering the 



coupling matrix expression (I10[) . For the p electrons and 
holes in the energy range 4 eV< e ^ 5 eV, both sets 
of data look very similar. Nevertheless, a more careful 
analysis (presented below) shows some differences which 
are reflected in the calculated inelastic lifetimes. 

The strongest SO effect in lifetime is observed for p 
states in vicinity of the BZ center. In Fig. [3l the arrow 
points to the r„ = 2,k results for k's vectors close to the 
r point calculated with the inclusion of the SO splitting 
in the band structure. For completeness, in Table [TTI the 
linewidth value for the lowest p band states at some k's 
along the TX direction is presented. These data demon- 
strate how the giant SO splitting of the p bands at the T 
point reduces the linewidth of the states in this band by 
as much as ~ 60%. This is mostly due to the reduction 
of phase space for final states. 
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Very recently lifetime measurements have been per- 
formed for the thin Pb films with (111) orientation grown 
on semiconductors . 41 ' 42 Such orientation corresponds to 
the quantization of bulk electronic states along the TL 
symmetry direction. Figure [4] presents the calculated 
damping rates for the excited electrons in p states along 
this direction. In difference with the situation for the 
electronic states around the T point one can see that the 
SO interaction produces small effect on lifetime for these 
states. On the other hand, the presented calculated re- 
sults are in good agreement with the experimental data 
for thin film a 41 ' 42 signalling that low energy quasipar- 
ticles dynamics in such systems can be well described 
by that in bulk Pb. Note also how the lifetime data for 
states along the TL direction in Fig. 0]both in the SC and 
SO calculations are well fitted by a quadratic function, 
r = a(E — Ep) 2 , over rather extended energy range. 



A. Dependence on the quasiparticle energy 

Although the damping rate T n k of a quasiparticle in a 
state (n, k) depends on its band index n and wave vector 
k, one can define T(e) as an average of r n ic over states 
with all wave vectors and bands lying at the same energy 
in the BZ. In Fig. [5] the results for r(e) are presented 
for quasiparticles in the p-like bands. The quasi-linear 
energy dependence of the averaged damping rate for holes 
with energies e 5j -2.5 eV requires a separate analysis of 
the dependence of r(e) for hot electrons and holes. 

In Fig. [5] different energy dependence of the averaged 
damping rate for holes is observed for energies below and 
above e s» -2.5 eV. At energies e above -2.5 eV, the aver- 
aged linewidth presents a quasi-quadratic dependence on 
quasiparticle energy. For more detailed analysis of this 




0.0 0.5 1.0 1.5 2.0 2.5 



Energy (eV) 

FIG. 4: (Color online) Calculated ab initio damping rates for 
excited electrons in states near the Fermi level with wave vec- 
tors along the TL direction excluding (squares, dashed line) 
and including (circles, solid line) SO interaction. Lines are 
fits by a quadratic function, F — a- e 2 . 



dependence, we fitted the T(e) data for -2.5 eV ^ e 
by a third-order polynomial: 

r fit (e) = a v e + a 2 - £ 2 + a r e 3 . (14) 

The coefficients Oj obtained from this fitting are pre- 
sented in Table HTT1 for all four kinds of linewidth calcula- 
tions. Note the highest value of the curvature in the SO 
case which is a direct consequence of the avoiding band- 
crossing effect produced by inclusion of the SO coupling 
in the Hamiltonian seen in Fig. Q] In the -4.2 eV ^ e ^ 
-2.5 eV energy range the coefficients ni w »2 » 03 for 
the SC, SO, LDF-SC calculated curves presenting strong 
contribution of a linear term with exception of the FEG 
results which are well described by a quadratic function 
for any energy. Hence the quasi-linear behavior of the 
averaged damping rates in all the calculations performed 
with the use of the ab initio eigenstates reflects the non- 
free-electron like behavior of holes in the lowest p band 
at energies below -2.5 eV. 

In the case of electrons, all curves in Fig. [5] present 
an apparent quadratic dependence on the quasiparticle 
energy. Nevertheless, in this case we performed also the 
fitting procedure with the use of expression (|14p. Ta- 
ble Mil presents the obtained corresponding coefficients 
a,i as well. At first sight, all four studied curves show the 
expected quadratic energy dependence with the FEG re- 
sults presenting the greatest deviation. At the same time, 
from the data of Table IIIII it is clear that the SO curve 
has the strongest curvature (bigger quadratic coefficient). 

An interesting point comes from the comparison of 
two T(e) curves calculated using the Lindhard dielec- 
tric screening (LDF-SC and FEG cases), with the two 




-8 -6 -4 -2 2 4 6 



Energy (eV) 

FIG. 5: (Color online) Energy dependence of averaged damp- 
ing rates, r(e), of quasiparticles in the p-like bands obtained 
in the scalar-relativistic calculation (dashed line), the one in- 
cluding the SO interaction (solid line) and the calculation 
which includes the ab initio eigenstates and the Lindhard di- 
electric function (dashed-dotted line). The dashed-dotted- 
dotted line shows the same as the orange solid line in Fig. 

El 
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TABLE III: Fitting coefficients of the T(e) data using a third-order polynomial (|14|l for holes in the energy interval above -2.5 
eV and p electrons. Meaning of abbreviations "SC", "SO", "LDF-SC", and "FEG" is explained in the text. a 2 in eV" 1 and 
as in eV -2 . 







holes 






electrons 






SC 


SO LDF-SC 


FEG 


SC 


SO LDF-SC 


FEG 


"1 

(12 

az 


-0.0030 
0.0195 
0.0000 


-0.0000 -0.0020 
0.0313 0.0150 
0.0050 -0.0009 


0.0000 
0.0216 
-0.0002 


0.0007 
0.0175 
0.0007 


-0.0020 0.0002 
0.0190 0.0180 
0.0004 -0.0007 


0.0060 
0.0161 
-0.0008 



TABLE IV: Effective charge density parameters obtained 
from Eq. I|15p on base of four sets of data as explained in 
the text. 5r a stands for the deviation from a conventional 
value for Pb r s = 2.298 a.u. 

SC SO LDF-SC FEG 



rf 2.37 2.49 2.18 2.33 

5rf(%) +3.0 +8.3 -5.2 +1.3 



TABLE V: Lifetime of excited electrons (holes) at four differ- 
ent values of |e|. All values are in fs. 



M 


SC 


SO 


LDF-SC 


FEG 


QF 


0.5 eV 


130(114) 


114(87) 


162(145) 


128(130) 


133 


1.0 eV 


34(29) 


31(23) 


41(36) 


34(31) 


33 


2.0 eV 


9(8) 


8(7) 


10(9) 


9(7) 


8 


3.0 eV 


4(3) 


4(3) 


4(3) 


5(3) 


4 



curves calculated using the ub initio screening (SC and 
SO cases). The former ones deviate considerably for ener- 
gies above ~ 3 eV from the two latter. Nevertheless, the 
curvature (i.e., 02 coefficients) is similar for all four curves 
(see Table llll)) . It is the negative sign of the 03 coefficient 
in the case of the averaged damping rates calculated with 
the Lindhard screening giving origin of the observed de- 
viation in the high-energy side. Hence, though a 2 ^03, 
for sufficiently high electron energies (e ^ ~ 3 eV) the 
cubic term in the dependence of T(e) on the energy can 
start to play significant role, the sign of 03 being related 
to the screening. 

From the results for T(e) an effective charge den- 
sity parameter rf can be derived with the use of the 
Quinn-Ferrel expression (fT3"|) . Fitting again the averaged 
linewidth curves with Eq. (|14p. new ai coefficients are 
calculated. Because of the approximations that lead to 
Eq. (flUf (see Refs. [H and l59l). Tq F is only valid for 
very low energy quasiparticles. Thus, the new fitting is 
carried in the energy range -1 eV ^ e ^ 1 eV, averag- 
ing the effect of possible different curvatures for electrons 
and holes. Finally, the rf parameters are found using 
the following expression: 

rf = (399.7 x a 2 ) 2/5 . (15) 

In Table IIVI the calculated rf and their deviation from 
the conventional r s = 2.298 a.u. for Pb are presented. 
Taking into account that fitting details, i.e., the choice 
of the polynomial or of the energy range applied in the 
fitting procedure, may influence the exact resulting val- 
ues of rf, a good agreement with r s is found. The value 
of the effective density parameter in the FEG case is the 
closest one whereas the higher value of rf for the SO 
case in comparison with the purely scalar-relativistic one 
means that the effective screening at such low quasipar- 
ticle energies is stronger in the former case. 



In Table [V] we present the calculated values for av- 
eraged inelastic lifetimes for excited electrons and holes 
in Pb at |e| = 0.5, 1.0, 2.0 and 3.0 eV obtained in all 
the calculations. For excited electrons, the QF, FEG, 
and SC calculations give similar results. Noting that the 
calculated lifetimes using the self-consistent eigenstates 
and the Lindhard screening differ from that three ones, 
we can conclude that, even if lead is a free-electron-like 
metal, lifetime of quasiparticles in bulk Pb is the result 
of a balance between screening and localization. Inter- 
estingly, the lifetimes in this energy interval calculated 
including the SO splitting in the band structure are the 
lowest ones for each quasiparticle energy, both for elec- 
trons and holes. 

As seen in Fig. [TJ at the W point only one unoccupied 
p band approaches but does not cross Ep, presenting a 
local minimum. This proximity to the Fermi level to- 
gether with the absence of available unoccupied states 
for decay with small momentum transfers leads to no- 
tably longer inelastic lifetimes for the states around the 
W point. This is a strong band structure effect. In Fig. [5] 
we compare the averaged inelastic lifetimes, r av , of elec- 
trons with very low energies (solid line) with the values 
calculated for states close to W (shown by symbols) ob- 
tained in the scalar-relativistic calculation (note that for 
this analysis the SO interaction has negligible effect). For 
comparison, the FEG results are shown by dashed line. 
As can be seen, upon approaching the W point (reducing 
energy) t increases faster in comparison with the aver- 
aged inelastic lifetime at those energies, being by s» 50% 
higher than r av at the same energy in other parts of the 
BZ. Hence, around the W point, r behaves in a non-free- 
electron-like manner. 
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0.4 0.5 0.6 

Energy (eV) 

FIG. 6: (Color online) Averaged inelastic lifetimes of elec- 
trons (solid line) and the lifetimes calculated for states at k's 
in vicinity of the W point (circles) evaluated in the scalar- 
relativistic calculation. Triangle highlights r at the W point. 
Dashed line: FEG results. 




1.0 1.5 2.0 2.5 3.0 3.5 

Binding energy (eV) 

FIG. 7: (Color online) Ratio r^jr^ as a function of the quasi- 
particle energy for four levels of calculations. The dashed line 
separates two energy regions with different physical behavior 
(see text). Lines are guides to the eye. 



B. Ratio of the lifetime of electrons and holes 

In Fig. [3 the ratio of lifetimes for electrons, r e , and 
holes, Th, at the same absolute value of the quasiparti- 
cle energy are shown for four different calculations. The 
curves show two energy ranges with different behavior. 
At energies |e| >2 eV all the four curves present a quasi- 
linear behavior. However, for quasiparticle energies be- 
low 2 eV, only the FEG curve remains quasi-linear. The 
other three calculations present quite different depen- 
dence of the r e /r/j ratio at those energies. The quasi- 
linear behavior of t g /t^ is found in the homogeneous elec- 
tron gas calculations (see ). Hence in bulk Pb the band 
structure effects are important in the electron-electron 



inelastic scattering processes for quasiparticles with en- 
ergies less than 2 eV. As the LDF-SC curve deduced from 
the lifetime results obtained with the use the Lindhard 
screening and the ab initio eigenstates presents also the 
band structure effects, these effects are the consequence 
of using the true eigenstates in the evaluation of the cou- 
pling matrices [see Eq. (10)], and not of the ab initio 
screening used. Note also that r e > 77,,V|e|, for all four 
levels of calculations. 



C. Effect of SO interaction on lifetime 

Modifications in the Pb band structure upon inclusion 
of the SO interaction increase lifetime of electrons in the 
d bands and reduce that of holes in the s ones, as ex- 
pected from Fig. [3J In all the cases the variation in r 
upon inclusion of the SO term does not exceed 10% in 
comparison with values obtained in the scalar-relativistic 
calculations. In Fig. [5] the ratio r(e )/t(£ ) as a func- 
tion of energy is plotted for the p holes and electrons. For 
holes with e — 2 eV and electrons with 2.5 eV ^ e 
5.5 eV inclusion of the SO splitting increases the quasi- 
particle lifetimes, in both cases less than 10%. On the 
other hand, for holes with binding energies less than 2 eV, 
spin-orbit induced splitting reduces lifetimes as much as 
by 30%. For electrons with e ^ 2.5 eV, as can be seen 
in Fig. [8j the SO interaction lowers r by w 10%. One 
can expect that this tendency will continue for electronic 
states around the Fermi surface. However, the k mesh 
used in the samplings over the BZ in present lifetime cal- 
culations does not allow to obtain reliable lifetime values 
for the states with binding energies less than 0.5 eV. 
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FIG. 8: (Color online) Ratio r(e acaiar )/r ( e so ) as a function 
of the quasiparticle energy for p-like states, both for holes 
(circles) and electrons (squares). The vertical line separates 
the quasiparticle energy ranges with different scales. 
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V. CONCLUSIONS 

We have presented GW-RPA calculations of the inelas- 
tic damping rates of low energy quasiparticles in bulk Pb 
and studied in detail the band structure effects as well as 
the consequences of the inclusion of the spin-orbit inter- 
action. A strong localization effect is found for electrons 
with energies 5.5 eV ^ e ^ 8 eV, where lifetime for d 
states is almost 2 times longer than that for p states. 

The states in the lowest valence p energy band at the 
center of the BZ reduce their damping rates by roughly 
~ 60% upon inclusion of the spin-orbital interaction in 
comparison with the scalar-relativistic calculation. For 
the p electrons the damping rates T(e) averaged over the 
BZ present a quadratic dependence on energy e, whereas 
for p holes the quasi-linear dependence of T(e) at e be- 
low -2.5 eV reflects a non- free-electron- like nature of the 
electronic states at the bottom of the lowest p band. The 



calculated lifetimes ratio r e (e)/Th(e) reveals that band 
structure effects are important in electron-electron inelas- 
tic scattering processes for quasiparticles with binding 
energies smaller than 2 eV. Present lifetime calculations 
for bulk Pb are in good agreement with inelastic damp- 
ing rates of quantum well states of metallic thin films 
measured very recently presenting thereby evidence that 
quasiparticle dynamics can be considered as being bulk- 
like even in very thin Pb films. 
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